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ABSTRACT 

We study the Spectral Energy Distributions, SEDs, (from FUV to MIR bands) of the first sizeable 
sample of 34 low-luminosity radio galaxies at high redshifts, selected in the COSMOS field. To model 
the SEDs we use two different template-fitting techniques: i) the Hyperz code that only considers 
single stellar templates and ii) our own developed technique 2SPD that also includes the contribution 
from a young stellar population and dust emission. The resulting photometric redshifts range from 
z~0.7 to 3 and are in substantial agreement with measurements from earlier work, but significantly 
more accurate. The SED of most objects is consistent with a dominant contribution from an old stellar 
population with an age ~ 1 — 3x 10 9 years. The inferred total stellar mass range is ^10 10 — I0 12 
Mq. Dust emission is needed to account for the 24 fim emission in 15 objects. Estimates of the dust 
luminosity yield values in the range Ldust ~ 10 43 5 — 10 45 5 erg s _1 . The global dust temperature, 
crudely estimated for the sources with a MIR excess, is ~ 300-850 K. A UV excess is often observed 
with a luminosity in the range ~ 10 42 - 10 44 erg s" 1 at 2000 A rest frame. 

Our results show that the hosts of these high-z low-luminosity radio sources are old massive galaxies, 
similarly to the local FRIs. However, the UV and MIR excesses indicate the possible significant 
contribution from star formation and/or nuclear activity in such bands, not seen in low-z FRIs. Our 
sources display a wide variety of properties: from possible quasars at the highest luminosities, to 
low-luminosity old galaxies. 

Subject headings: Galaxies: active - galaxies: high-redshift - Galaxies: nuclei - Galaxies: evolution - 
Galaxies: photometry - Infrared: galaxies 



1. INTRODUCTION 

The search for radio-loud Active Galactic Nuclei 
(AGN) is one of the widely used tools to study the dis- 
tant universe (z > 1). Indeed the fir st objects found 
at z > 1 were radio-loud quasars (see iStern fc Spinradl 
119991 and references therein) and a large number of radio 
galaxies (~300) are know n to exist at high redshifts (e.g., 
IMilev fc De Breuck|[2008h . 

Flux-limited samples of radio galaxies such as the 3CR 
and its deeper successors 6C and 7C catalogs are affected 
by a tight redshift-luminosity correlation. This, along- 
side with the steep luminosity function of radio sources, 
gives rise to a selection bias which results in the pres- 
ence of high power sources predominantly at high red- 
shifts only and low power sources exclusively at low red- 
shift. Therefore, our knowledge of radio galaxies and 
their hosts at high redshifts is ex clusively based on stud- 
ies of powerful (edge-brightened, iFanaroff k, Rilevil 19741) 
FR lis. Thus, the properties of the low-luminosity AGN 
population in the early Universe is essentially unknown. 
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The missing pieces of the puzzle might be obtained by 
studying a sample of low luminosity ('edge-darkened', 
IFanaroff fc Rilevl 11974ft FR Is at high redshifts. The 
lower power of these objects with respect to that of high- 
z FR lis might allow us to analyze with more clarity the 
properties of the hosts at high redshifts. Furthermore, 
in terms of host galaxy properties, low power sources are 
both more abundant and most likely more similar to qui- 
escent galaxies (in terms of host galaxy properties) than 
FR lis. 

Only a few F R Is at high redshift are present in 
the 7C sample dHevwood et al.l 12007ft and two wer e 
possibly found i n HD F North (jSnellen fc Best! 12001ft . 
iChiaberge et al.l ()2009h obtained the first sizeable sam- 
ple of 37 FR I candidates at z > 1, located in the 2 
deg 2 area of the sk y observed by the COSMOS survey 
(|Scoville et al.|[2007|) . In this paper we perform a detailed 
analysis of the photometric properties of the host galax- 
ies of these radio sources, by taking advantage of the 
large multi-wavelength coverage provided by the COS- 
MOS collaboration. This allows us to derive their spec- 
tral energy distributions (SEDs) from ~ 0.15 /im to 24 
(Ltm. We then model the obtained SEDs with synthetic 
stellar population templates with the aim of both infer- 
ring the properties of each galaxy and of deriving an 
a ccurate estimate o f the photo metric redshi f ts. 

iMobasher et~aT1 ()2007ft and lllbert et al.l (|2009ft per- 
formed a similar study with the aim of determining the 
photo-z for the whole COSMOS samp le. However, as 
already noted by IChiaberge et al.l (|2009ft , the faint coun- 
terparts of these radiogalaxies might be easily misidenti- 
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fied or be missed from the COSMOS catalog. Therefore 
a careful object-by-object study of the sample is nec- 
essary to correctly identify the genuine counterparts of 
the radio sources at all wavelengths. This allows us to 
obtain reliable SEDs and perform a robust study of the 
hosts of these objects. In this work, using the multi-band 
data provided by the COSMOS survey (described in Sec- 
tion[3]), we obtain the photometric data from the catalog 
checking visually the correct identification of each object 
of the sample (Section 0]). In Section [5] we describe the 
codes used to model the SEDs, Hyperz and 2SPD. We 
present the results obtained from the SED modeling in 
Section|BJ the photometric redshifts, the radio power dis- 
tribution, and the host properties, such as stellar ages, 
masses and dust and UV components. In Section [7] we 
summarize the results and we discuss our preliminary 
findings. 

We adopt a Hubble constant of Ho = 71 km s _1 Mpc -1 , 
fi TO = 0.27 a nd VL vac = 0.73, as given by the WMAP cos- 
mology (e.g., ^E^^A^QQMMo^^^MiM- A11 
the magnitudes are in AB mag system, if not otherwise 

specified. 

2. SAMPLE 

iChiaberge et al.1 (|2009| ) selected a sample of 37 high-z 
low-power radio galaxies at z ~ 1-2 in the Cosmic Evolu- 
tion Survey (COSMOS) field (IScoville et al.ll200l . This 
is the first sizeable sample of low-power radio sources 
at mid- to high-redshifts. This sample has been ob- 
tained with a four-steps multi-waveleng th selection pro- 
cess, which is described in detail in IChiaberge et al.l 
(|2009f ). Here we briefly summarize the selection pro- 
cedure which depends on the following assumptions: i) 
the FR I/FR II break luminosity at 1.4 GHz does not 
change with rcdshift, and ii) the host properties of dis- 
tant low-luminosity sources are similar to those of high- 
power FR lis in the same redshift range. 

The first step con sists of selecting ra dio sources from 
the FIRST catalog (jBecker et al.l fl995h whose 1.4 GHz 
total flux corresponds to the luminosities expected by 
FR Is at 1< z <2. The second step is based on a ra- 
dio morphological classification: the radio sources which 
show clear "edge-brightened" structures are rejected in 
order to exclude the bona-fide FR lis from the sample. 
The third step implies the optical identification of the ra- 
dio sources in the COSMOS optical images. The objects 
associated with host galaxies brighter than i < 22 were 
rejected to exclude lower redshifts starburst galaxies. As 
a final step, u-band dropouts are also rejected since these 
are objects most likely located at redshift higher than 
z=2. The resulting sample consists of 37 FR I candi- 
dates. A posteriori, the photometric r edshifts range of 
most of them is betwe en ~1 and ^2 ijMobasher et al.l 
120071: lllbert et al.ll2009f) . with the exception of 3 objects 
(namely, 7, 27, and 66 7 ) that we exclude from any further 
analysis because out of the redshift range of interest. We 
are then left with a final sample of 34 objects. The list 
of the radio sources is presented in Table [TJ 

3. DATA 

3.1. UV, optical, and IR photometric data 

7 The 'naming convention' of the sources used in the paper is 
based on the listing number in the COSMOS catalog. 



The photometric data used to build the SEDs of 
our sources are ta ken from the COSMOS survey 
(|Scoville et al.l 12001 . The survey comprises ground 
based as well as imaging and spectroscopic observations 
from radio to X-rays wavelengths, covering a 2 deg 2 area. 
Given the high sensitivity and resolution of these data, 
COSMOS provides samples of high redshift objects with 
greatly reduced cosmic variance as compared to earlier 
surveys. 

Ground-based UV, optical, and IR observati o ns an d 
data reduction ar e pr esented in Capa k~et al.l (|2007| ). 
iCapak et all (j2008l ) and iTaniguchi et all (|2007K A mul- 
tiwavelength photometric cata l og w as generated using 
SExtractor (|Bertin fc Arnoutsl fl996h . The COSMOS 
catalog is derived from a combination of the CFHT i* 
and Subaru i + images, to which the authors refer as 
'I-band images. The catalog includes objects with to- 
tal ("mag-auto") I < 25 and searches for counterparts 
in a radius of l"around the I-band detection. At fainter 
magnitudes the catalog begins to be incomplete and have 
more spurious detections, and photometric redshifts are 
poorly constrained. 

For this study we use the COSMOS Intermediate 
and B road Band Photometry Catalog 2008 (ICapak et al.l 
|2008| ) 8 which provides the multiwavelength magnitudes 
of our sources from FUV to K bands. Narrow band fil- 
ters are not considered due to the possible strong con- 
tamination of emission lines from the AGN and their low 
signal-t o-noise ratio. 

HST dKoekemoer et alj|2007l) and S pitzer data (both 
IRAC and MIPS. iSanders et aLll2007| ) are also included 
in the COSMOS survey. The later is presented in two 
separate catalogs. S-COSMOS IRAC 4-channel Photom- 
etry Catalog June 2007 is for IRAC data. S-COSMOS 
MIPS 24 Photometry Catalog October 2008 and S- 
COSMOS MIPS 24 /im DEEP Photometry Catalog June 
2007 are for MIPS data at 24 fj,m with two different flux 
limits, 0.15 and 0.08 mJy, respectively. 

For the sake of clarity, we summarize here the available 
datascts from the COSMOS survey used in this work 
(Table ©: 

1. GALEX: UV dat a were taken using GALEX 
(jMartin et al.ll2005T ). The observations, performed 
as part of the GALEX Deep Imaging Survey, are 
in the NUV and FUV bands, respectively with an 
angular resolution of 5"6 (NUV) and 4'.'2 (FUV) 
(jMorrissev et al.ll2007D . 

2. HST: the COSMOS HST data are single-orbit 
F814W ACS images. Th ey have the highest an- 
gular resolution (~ 0"09. iKoekemoer et al.l 120071 ) 
among the COSMOS images. 

3. SUBARU: The Suprime-Cam instrument on the 
Subaru telescope observed the COSMOS field in 
six broad bands (Bj, g + , Vj, r + , i + , z + ) with a n 
angular resolution of ~0 / /2 (jTaniguchi et al.ll2007t l. 

4. CFHT: the Canada-France-Hawaii Telescope 
(CFHT) provides u* and i* images using Megacam 
jBoulade et al.|[2003h and K images with Wircam. 

8 http:/ /irsa. ipac.caltech.edu/cgi-bin/Gator/nph-dd 
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TABLE 1 

Radio positions, redshifts and radio flux of the sample 



ID 


RA 




DEC 


z phot,Ilbert 


Zspec 


Fnvss 


Ffirst 


1 


150.20744 


2 


2818749 


QO +um 
u - yz -0.06 


0.8827 a -0.8823 b 


<2.5 


1.79 


2 


150.46751 


2 


7598829 




2.6 


1.08 


3 


150.00253 


2 


2586310 


J " MD -0.41 




5.2 


4.21 


4 


149.99153 


2 


3027799 


1.45«;f 5 




7.5 


5.99 


5 


150.10612 


2 


0144780 


1.84«-f 2 




3.4 


1.30 


11 


150.07816 


1 


8985500 


1 01 +0.60 
— 24 




<2.5 


1.13 


13 


149.97784 


2 


5042069 


1 09+° : ° 6 

— 06 




2.4 


1.51 


16 


150.53772 


2 


2673550 


n qk+0.05 
U - MO -0.02 


0.9687 6 


4.4 


5.70 


18 


149.69325 


2 


2674670 


0.92+°-?,i 

^ ^ — 0.09 




5.1 


4.39 


20 


149 £'3909 


2 








<2.5 


1.33 


22 


149 89^08 


2 


fi9991 44 




<2.5 


2.74 


25 


150.45673 


2 


5597000 


^ -0.03 


0.7917 6 * 


2.7 


2.18 


26 


149.62114 


2 


0919881 


x ' iu -0.02 




3.2 


1.88 


28 


149.60064 


2 


0918673 




2.4 


1.77 


29 


149.64587 


1 


9529760 


1 r.n+0.45 

— 30 




2.3 


2.12 


30 


149.61542 


1 


9910541 






2.4 


1.26 


31 


149.61916 


1 


9163600 




0.9132 a -0.9123 i ' 


4.1 


3.71 


32 


149.66830 


1 


8379777 




3.1 


1.31 


34 


150.56023 


2 


5861051 


1 4 n+0.65 
^ -0.32 




4.5 


5.25 


36 


150.55662 


1 


7913361 




3.3 


3.19 


37 


150.74336 


2 


1705379 


-r nrr + 09 




2.2 


1.87 


38 


150.53645 


2 


6842549 


-, -. r» + 10 

!■ 12 0.05 




11.6 


10.01 


39 


149.95804 


2 


8288901 


, n8 +0.03 




<2.5 


1.37 


52 


149.90590 


2 


3964710 


n 7 4+0.02 
u ' ' 4 -0.03 


0.7417 6 


<2.5 


1.54 


70 


150.61987 


2 


2894360 


9 91 +0.57 
z ' z -0.37 




4.5 


3.90 


202 


149.99506 


1 


6324950 


1 iq+°- 24 
i - i9 -0.14 




3.3 


1.08 


219 


150.06444 


2 


8754051 






<2.5 


1.85 


224 


150.28999 


1 


5408180 


i - iu -0.03 




3.2 


3.31 


226 


150.43864 


1 


5934480 


x - '°-0.14 




<2.5 


1.19 


228 


149.49455 


2 


5052481 


1 8Q+ - 74 
i - 89 -0.59 




3.7 


2.04 


234 


150.78925 


2 


4539680 


i io+ - 12 

1 - iu -0.03 




5.2 


4.43 


236 


150.70554 


2 


6296339 


2.132 c 


7.0 


7.10 


258 


149.55934 


1 


6310670 


90+ 02 
u ' yu -0.02 


0.9009 a 


3.7 


2.24 


285 


150.72131 


1 


5823840 


-, 91 +0.06 
± ' zl -0.08 




3.5 


2.95 



Note. — Colu mn description: (1) ID number of the object; (2)-(3) right ascension and d eclination of rad io source; (4) photometric 
redshifts foun d by Ilbert et al. 12009 ); (5) spect r oscopi c redshift: a from zCOSMOS survey (Lillv ct al. 2007), 6 form Magellan survey 
I'lhimp ct al. 20 03), and c from [Frescott ct al. (200TJ). The spectroscopic redshift of object 25 (marked with *) is considered incor- 
rect (see Section ISTU ; (6)-(7) flux at 1.4 GHz from NVSS (from http://www.cv.nrao.edu/nvss/NVSSlist.shtml) and FIRST (taken from 
Chiabcrge ct al. 2009) catalogs in mjy. 



5. UKIRT: near infrared Wide Field camera (WF- 
CAM) on the United Kingdom Infrar ed Telescope 
(UKI RT) provides the J-band images (|Casali et al.l 
l200l . 

6. NOAO: Ks data are taken at the Kitt Peak 
National Observatory (KPNO) telescope with 
FLAMINGOS and The Ccrro Tololo Inter- 
American Observ atory (CTIO) telescope 
(jCapak et al.l l2007f l. These telescopes belong 
to the National Optical Astronomy Observatory 
(NOAO). 

7. Spitzer: Spitzer cycle-2 S-COSMOS is an in- 
frared imaging sur vey of the COSMOS field 
([Sanders et al.ll2007h . They obtained observations 
with the IRAC camera in four channels, at 3.6, 4.5, 
5.6, and 8 fxm, and with MIPS in the 24, 70, 160 /xm 
band. We only consider the data at 24 /zm, since 
no object is detected at longer wavelengths (with 
the exception of object 37 which is also detected at 



70 /im). 

3.2. Radio data 

IChiaberge et al.1 (|2009ft selected the sample in the 
COSMOS field using radio fluxes at 1.4 GHz from the 
FIRST survey (jBecker et al.|[l995T) . The data are ob- 
tained with the VLA in B configuration with an angular 
resolution o f ~5 ,r and reach a flux limit of ~1 mjy and 
are listed in IChiaberge et al.l (|2009| ). 

In addition, the NVSS survey (|Condon et al.l 119981) 
provides 1.4 GHz radio data for our sample, but at lower 
resolution (~45 ") and with a higher flux density limit 
(~2.5 mjy) than the FIRST survey. These differences 
imply that seven of our objects are missed in the NVSS 
catalog. Nonetheless, the NVSS data arc useful since 
they are more sensitive to diffuse low surface brightness 
radio emission that the FIRST data. The radio fluxes are 
taken from the NVSS archive 9 . A NVSS/FIRST compar- 

9 http:/ /www. cv.nrao.edu/nvss/NVSSlist.shtml 
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TABLE 2 

COSMOS BROAD BANDS AND THEIR PROPERTIES. 



Filter 


Telescope 




FWHM 


sensitivity 


FUV 


GALEX 


1538.6A 


230.8A 


25.7 


NUV 


GALEX 


2315. 7A 


789.1 A 


26.0 


u* 


CFHT 


3911.0A 


538.0A 


26.5 


Bj 


Subaru 


4439.6A 


806.7A 


27.0 


9+ 


Subaru 


4728.3A 


1162.9A 


27.0 


Vj 


Subaru 


5448.9A 


934.8A 


26.6 


r+ 


Subaru 


6231.8A 


1348.8A 


26.8 


i* 


CFHT 


7628.9A 


1460.0A 


24.0 


i+ 


Subaru 


7629.1 A 


1489.4A 


26.2 


F814VF 


HST 


8037.2A 


1539.0A 


27.2 


z+ 


Subaru 


9021.6A 


9021.6A 


25.2 


.1 


UKIRT 


12444.1 A 


1558.0A 


23.7 


K S 


NOAO 


21434.8A 


3115. OA 


21.6 


K 


CFHT 


21480.2A 


3250.0A 


23.7 


IRAC1 


Spitzcr 


35262. 5A 


7412.0A 


23.9 


IRAC2 


Spitzcr 


44606.7A 


10113.0A 


23.3 


IRAC3 


Spitzcr 


56764.4A 


13499.0A 


21.3 


IRAC4 


Spitzcr 


77030.1 A 


28397.0A 


21.0 


MIPS1 


Spitzer 


23.68/xm 


4.7/^m 


29.6 



Note. — List of the filters used for our analysis, the associated 
telescopes, the effective wavelengths A e jj, the PSF full- width half 
maximum (FWHM) of the images in each band, and their sensi- 
tivities at 5<t (mag). 

ison indicates that the flux ratio between the two catalogs 
is usually between <~ 1 and 2, and never larger than ~ 3. 
Table Q] shows the NVSS and FIRST radio fluxes of the 
objects. 

3.3. Spectroscopy 

The COSMOS survey provides spectroscopic data 
from the Very Large Telescope (VLT) (zCOSMOS, 
iLillv et all 120071) and fro m the Magellan (Baade) tele- 
scope (jTrump et al.ll2007t ) 

zCOSMOS is a large-area redshift survey which con- 
sists of two parts. The first, namely zCOSMOS-bright, 
considers a magnitude-limited (I-band mag <22.5) sam- 
ple of ~20,000 galaxies located within the central 1 deg 2 . 
Spectra cover the wavelength range 5500 A < A < 9000 
A. The second, namely zCOSMOS-deep, is an ongoing 
survey (not yet public) of ^10,000 blue galaxies in the 
same field filtered with a color selection to be in the range 
of 1.4 < z < 3.0. The spectra cover the wavelength range 
of 3600 A < A < 6800 A. 

The Magellan survey presents spectroscopic redshifts 
for the first 466 X-ray and radio-selected AGN candi- 
dates. The wavelength coverage of these spectra is ~ 
5500-9200A. Their redshift yield is 72% for i AB < 24 
and >90% for i AB <22. 

In this work we use only spectroscopic measurements 
with a confidence level greater than 99%. Six objects of 
the 34 objects are included in the spectroscopic surveys 
described above with the required quality (Table [1]) . 

4. MULTI-BAND COUNTERPARTS IDENTIFICATION 

The process of counterparts identification of our ra- 
dio galaxies in the COSMOS catalog suffers from the 
limitations typical of any multiband survey (such as 
misidentification of targets with close neighbor or the 
contamination by nearby bright sources). We then pre- 
fer to perform our multi-band counterpart identification 



on each source by visually inspecting its multiband im- 
ages, rather than blindly use the data provided by the 
COSMOS catalog. 

We start the process looking for a I-band counterpart 
to the radio source in the COSMOS catalog by adopt- 
ing 0'.'3 as search radius. 29 identifications are found in 
our sample, most at distances smaller than O'.'l. Three 
objects (22, 28, and 32) are clearly visible in the I-band 
images, but they are not found in the COSMOS catalog, 
since they are below its detection threshold. In addition, 
there are two objects (2 and 36) for which no I-band 
counter part is found. 

For the I-band detected sources, we search in the COS- 
MOS broadband catalog that provides photometry (from 
the FUV to the K band) over an aperture of 3" diam- 
eter. The COSMOS catalog associates its counterparts 
in the remaining bands with the brightest and closest 
sources to its I-band position within a radius of 1". We 
extend this method including also the Spitzer/IRAC and 
MIPS catalogs, by using a larger search radius (2") due 
to their coarser resolutions with respect to the COSMOS 
broadband images. For the 5 sources not present in the 
COSMOS catalog (because they are too faint in the I- 
band) we note that in all cases a clear counterpart is 
instead present at longer wavelength providing us with a 
robust multiband identification. 

During the visual check of the individual counterparts, 
two main problems emerge in the COSMOS broadband 
catalog: i) the presence of nearby sources, within the 3" 
radius used for the aperture photometry, contaminating 
the broad band measurements of the genuine emission 
from the radio-galaxies of our sample; ii) the counter- 
parts to the i-band object does not always correspond to 
the same object over the various bands. Object 28 is a 
clear example of this situation (see Fig. [T]): from the J 
band to 8/im the brightest source in the field is coincident 
with the radio source, while in the i + and z + bands this 
is out-shined by a nearby object, causing an erroneous 
identification in the COSMOS catalog across the various 
bands. Furthermore, the presence of this neighbor also 
causes a strong (and dominant short-ward of the J band) 
contamination to the genuine flux of the radio-galaxy. 
Other examples of contamination and mis-identification 
are reported in FigO 

In order to amend these problems, we first perform a 
new 3" aperture photometry properly centered on the 
position of the radio source, to isolate the genuine emis- 
sion of the counterparts. In case of contamination from 
nearby source(s), we subtract from the flux resulting 
from the photometry centered on the radio source the 
emission from the ncighbor(s), limiting to the fraction 
that falls into the 3" aperture. When we cannot sepa- 
rate a source from a close source or when the counterpart 
is not visible in a given band, we measure an 1-er upper 
limit to the flux. 

At 24 /jm, when the counterpart of a source is not seen 
at the radio position, we use the detection limit of MIPS 
catalogs as upper-limit to the source flux. However, when 
a nearby source contaminates the emission of our target, 
we prefer to measure the upper-limit directly on the im- 
age at the radio position. In addition, in the NUV/FUV 
band, if the source is not detected in GALEX we prefer 
not to include upper limits in our analysis, because its 
corresponding flux is substantially higher than those at 
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Fig. 1. — Images of object 28 taken from the COSMOS survey. From upper left to lower right panels we show the images from Subaru 
i + , Subaru z+, UKIRT J, CFHT K and the four channels of Spitzcr/IRAC. The size of the image is 5'.'5 X 5'.'5. We overlap the radio 
contours (green lines) at 1.4 GHz from VLA COSMOS data on the frames to show the position of the radio source. 



larger wavelengths. 

In some cases the GALEX photometry returns appar- 
ently incorrect results. In Fig. [3] we show as an example 
the HST and GALEX images of object 13; this source 
is labeled as detected by the COSMOS catalog in both 
the FUV and NUV GALEX bands, while clearly there is 
no emission above the background level. The opposite is 
seen in MIPS measurements, when clearly visible objects 
are missing from the catalog (see Fig.|4|). In these cases, 
we do not consider the GALEX data, while for MIPS we 
provide a new estimate to the 24 /im flux. 

We also find an apparent error in the reported UKIRT 
J-band photometric points. We obtain our own measure- 
ments on a sample o f stars present in the COSMOS field 
(|Wright et al.l 12010ft using the PSF-matched and PSF- 
original J-band images 10 and we compare them to the 
magnitudes given by the catalog. This test reveals that 
the zero-point mag of the PSF-matched J-band images is 
higher than the corresponding value of the PSF-original 
J-band images, taking into account of the effect of the 
PSF convolution. This systematic difference of the zero- 
point magnitudes is AJ = 0.90±0.06, intriguingly simi- 
lar to the offset between Vega and AB mag system for 
the UKIRT J-band (3 A b = Jvcga + 0.94 for UKIRT, 
IHewett et al.l |2006f ) . We obviate this problem by per- 
forming 3"-aperture photometry on the J-band counter- 
parts on the UKIRT images smoothed with a Gaussian 
kernel to a 1.2-1.5" FWHM to reproduce the effect of the 
PSF matching. 

Once we obtain the correct aperture photometry of the 
various multi-band counterparts of the radio galaxies, we 
apply the appropriate aperture corrections. While this 

10 The PSF-matched images are obtained by convolving each 
PSF-homogcnized image with a Gaussian kernel that produced the 
same flux ratio between a 3"and 10" aperture for a point source to 
avoid PSF-matching problems in the multiband catalog. The PSF- 
original images are the pure images without any PSF matching. 
Both the images are provided by the COSMOS catalog. 



is not needed for the Subaru, CFHT, and UKIRT, for 
GALEX data we apply an aperture correc tion by mul- 
tiplyi ng the total flux by a factor of 0.759 (jCapak et al.l 
[2007h . Similarly, since we select from the IRAC catalog 
the 2'.'9 aperture photometry (because it is closest to the 
3" aperture of the COSMOS broadband catalog among 
the different optional apertures) , the IRAC aperture flux 
was converted to total source flux by using the aperture 
correction factors (from the IRAC manual) of 1.19, 1.27, 
1.48, and 1.27 for the four channels. The MIPS data in- 
stead already include the aperture correction. Since the 
CFHT telescope is more sensitive and with higher reso- 
lution than the images from the NOAO telescopes, this 
usually results in far smaller errors for the K-band mag- 
nitudes. In such cases, we prefer to use only the CFHT 
K-band data. The results of the corrected 3"-apcrturc 
photometric measurements arc tabulated in Tables [5] and 

m 

5. SED FITTING 

The SEDs are derived by collecting multiband data 
from the FUV to the MIR. Since not all of the objects are 
detected in the entire set of available bands, the number 
of detections used to constrain the SED fitting ranges 
from 15 to 19. 

The synthetic stella r templates used to model the ob- 
served SEDs are the iBruzual fc Charlotl pOOl ) (priv. 
comm.) and iMarastonl ( 20051 ) templates. These tem- 
plates are defined with different Initial Mass Fun ction 
(IMF) (|SalDeterlll955l IKrounall2l¥)ll : IChabrierll2003D with 
solar metallicity (Z Q = 0.02). These libraries contain 
composite stellar population (CSP) computed with dif- 
ferent star formation histories: a constant star-forming 
system (with constant star formation rate of 1 M Q /yr); 
a single burst of star formation; and ten /i-models 11 with 

These templates correspond to synthetic spectra computed 
with exponentially decaying star formation rate, if>{t) <x e' — */ r ) 
where r is the star formation timcscale. 
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Fig. 2. — Left panel: Subaru image of object 34 with superposed the radio contours (green) and the 3" aperture centered on the 
radio source (blue). This is a case in which, in addition to the radio-galaxy, other objects lie inside the integration aperture of the target, 
contaminating the COSMOS catalog flux measurement. Middle and right panels: Subaru z + and CFHT u* images of object 5 with 
superposed the radio contours. This is a case of misidentification of a source with a neighbor included within the 1" from the target) in 
different bands. The size of the images is 2"x 2". 



HST/ACS F814W GALEXNUV GALEX FUV 



HST/ACS F814W 



MIPS 24 micron 



Fig. 3. — From left to right: HST/ACS F814W, GALEX NUV and FUV bands images of object 13. The radio source (whose position is 
marked by the circles at the images center) does not apparently correspond to any UV emission, contrarily to the detection found in the 
catalog. The size of the images is 7'.' 5x7'.' 5. 

giant branch (TP-AGB) phase. These stars significantly 
contribute to the IR emission for ages higher than ~1 
Gyr. The resulting effect is that the Bruzual & Char- 
iot templates ar e bluer and give slightly older ages than 
iMarastonl (|2005f) models. 

5.1. SED fitting method: Hyperz 

We firstly estimate the photometric redshift for our 
sample by using the template fitting technique, Hy- 
perz (Version 1.3 , obtained fro m M. Bolzonella i n priv . 
comm.), described in detail by [Bolzonel la et al.l (|2000() . 
This SED fitting procedure is based on reproducing the 
overall shape of the SEDs and recognizing strong spectral 
properties, such as the 4000 A Ca break or Lyman break 
at 912 A. This code consists of a convolution of tem- 
plates, which represent the rest-frame SEDs for galaxies 
with different star formation histories, with the filter re- 
sponse functions given by the COSMOS survey. 

The basic assumption of this SED modeling is that 
the host galaxies of the radio sources in our sample are 
dominated by stellar emission. Thus the MIPS 24 fj,m 
photometric point is excluded because this emission is 
obviously not of stellar origin. Furthermore, object 236, 
a spectroscopically-confirmcd QSO (see Sect. I5.3[) . is 
excluded from this analysis. The code also takes into 



* 










Fig. 4.— Images of object 70 from ACS (left) and MIPS 24,1ml 
(right). This provides an example of a source clearly emitting 
at 24/xm, but not present in the COSMOS Spitzer/MIPS catalog. 
The position of the radio source is represented by the green radio 
contours from VLA COSMOS. The size of the images is 28"x28". 

time-decays of 0.1, 0.2, 0.3, 0.6, 1.0, 2.0, 3.0, 5.0, 10.0, 
and 15.0 Gyr. These templates include 221 tracks of ages 
from 0.1 Myr to 20 Gyr and cover the wavelength range 
from 91 A to 160 /im. Fig. [3] shows the stellar templates 
used fo r SED modeling. 

The iBruzual fc Chariot! 1)20091) models differs from 
those of IMarastonl ( 20*05f) mainly because of a different 
recipe for the inclusion of thermally pulsing asymptotic 
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Log ( A A ) Log ( K A ) 

Fig. 5.— Examples of SED templates from BC09 JBruzual fc CharloHl20fM) (left panel) and MA05 l lMarastonll2005l ) (right panel) in units 
of Lq A -1 . The templates correspond to different ages and different star formation histories. Arranging the BC09 and MA05 templates 
with increasing luminosity at 3550 A (following the dashed line), the different plotted models correspond to £t=15 Gyr decaying CSP of l.f 
Myr (red line); single stellar population of 20 Gyr (gray line); fi=0.l Gyr decaying CSP of fO Gyr (dark green line); fi=l0 Gyr decaying 
CSP of 2.8 Myr (light green line); (j,=0.3 Gyr decaying CSP 5.25 Gyr (orange line); fi=5 Gyr decaying CSP of 6.9 Myr (blue line); fj,=3 
Gyr decaying CSP of 17 Myr (yellow line); fi=0.6 Gyr decaying CSP of 1.7 Gyr (pink line); fi=2 Gyr decaying CSP of 38 Myr (fuchsia 
line); and p.=l Gyr decaying CSP of 0.2 Gyr (light blue line). We invite the reader to download the color version of the plot to distinguish 
the different lines. 



account extinction (by using the reddening law from 
iCalzetti et al.ll2000D which is applied to the templates. 
The reddened templates are then shifted in wavelength 
searching for the correct redshift. The best fit is obtained 
through x 2 minimization. The grid spacings in redshift 
and in Ay are 0.01 and 0.1, respectively. The Hyperz 
code gives the probability P(Xnd) associated with Xrid 
(i.e., the reduced x 2 7 which is defined as X 2 n i n l v ^ where v 
is the number of degrees of freedom). Because we are es- 
sentially interested in the redshift (the other parameters 
are affected by the degeneracy), the degrees of freedom 

are N/ i;ters - 1 

Similarly to lllbert et al.l (|2009f) , we increase the flux 
errors by 10% when the probability associated with the 
is less than 0.01. This factor does not shift the 



best-fit photo-z value but broadens the x 2 peak and de- 
rived redshift uncertainty The program provides the 
redshift of the object measured at the different confi- 
dence intervals denned by the values of xlid- We con- 
servatively choose the s olution at a 99% confidence level. 
iBolzonella et al.l <|2000f ) discussed the effects of degener- 
acy in the parameter space defined by star formation 
history, age, and reddening. However, they show that 
this does not strongly affect the value of the photometric 
redshift. Hyperzmass, a code which works in a similar 
way to Hyperz, returns the mass of the stellar popula- 
tion corresponding to the best fit and we adopt again the 
99%-confidence level to quote the errors on the mass. 

Fig. [15] (left panels) shows the SED fit obtained with 
Hyperz. The resulting properties from the modeling are 
listed in Table [3] Inspection of this figure indicates that 
Hyperz does not always return satisfactory results; in- 
deed, in several cases, it cannot reproduce the bluest part 
of the spectrum (possibly due to a young stellar com- 
ponent) or the substantial SED "bump" at long wave- 
lengths, clearly suggestive of dust emission and often ex- 
tending even to the IRAC 3.6 (im measurements. 

5.2. '2SPD' fitting technique 



In order to improve the quality of the SED fitting we 
developed a code that includes 2 Stellar Populations and 
Dust component (2SPD). More precisely, we take into 
account two different stellar populations, typically, one 
younger and one older (YSP and OSP, respectively). We 
model the dust component with a single (or, in a few 
cases, two) temperature black-body emission. Thus, at 
this stage, we can include the 24 /jm MIPS flux in the 
fitting process. 

The synthetic models used are the same as described 
in the previous section, but limiting to those with sin- 
gle stellar population with ages ranging from 1 Myr to 
12.5 Gyr. We adopt a dust-screen model for the extinc- 
tion normalized wit h the free parameter Ay, and the 
ICalzetti et all (|2000l ) law. 

The code searches for the best match between the sum 
of the different components (young and old stellar tem- 
plates and dust emission) and the photometric points 
minimizing the appropriate x 2 function. 2SPD returns 
the following free parameters: z, Ay, the age of the two 
stellar populations, the temperature of the dust com- 
ponents), and the various normalization factors. From 
these quantities we measure the stellar mass content of 
the two stellar populations at 4800 A rest frame. Sim- 
ilarly to the case of Hyperz, caution should be exerted 
before associating these values to physical quantities be- 
cause of degeneracy in the parameter space, apart from 
the photometric redshifts. 

To estimate the errors on the photo-z and mass deriva- 
tions, we measure the 99%- confident solutions for these 
quantities. This is computed by varying the value the pa- 
rameter of interest until the x 2 value increases by A% 2 = 
6.63, corresponding to a confidence level of 99% for that 
parameter. 

Fig. [15] (righ panels) shows the outputs of 2SPD code, 
while Table U presents the resulting parameters of the fit. 

The dust emission is usually poorly constrained. In 
many cases no excess above the stellar component is re- 
quired at A < 24/im. In this case dust is needed only 
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to account for the 24 fim data-point and often this mea- 
surement is an upper limit. The values reported in Ta- 
ble 2] correspond to the best fitting model for each galaxy 
and must be interpreted as approximate rather than real 
measurements of the dust emission properties. We will 
return to the dust properties in more detail in Sect. 16.41 

5.3. COSMOS photometric redshifts 

As a comparison with our results we also collect the 
photo metric redshifts f rom t he COSMOS rcdshift cat- 
alog. Mobash er et al.l (|2007l ) provide photometric red- 
shifts for 860,000 galaxies in the COSMOS field with i+ 
< 25. Their technique is based on a \ 2 template-fitting 
procedure applied to the SEDs derived with up to 16 pho- 
tometric points from the u to the K band. To evaluate 
the reliability of the derived photometric redshifts, they 
compare them with spectroscopic redshifts from a sam- 
ple of 868 galaxies in zCOSMOS with z < 1.2. The rms 
scatter between photometric and spectroscopic redshifts 
is o-(z P hot-zspec)/(i+ zsp ec) = 0.031 with a small fraction 
o f outliers (2.5%). 

Illbert et all (120091 ) substantially improve this analysis 
by including data from 30 broad, intermediate, and nar- 
row bands covering the wavelength range from UV to 
MIR. Redshifts are computed for 607,617 sources (with 
i + < 26). The method used by these authors accounts 
for the contribution of emission lines to the SEDs. Com- 
parison with 4,148 zCOSMOS spectroscopic redshifts in- 
dicates a dispersion of <T( Z phot-zspec)/(i+zs P ec) = 0.007 at 
i\ B < 22.5. Nevertheless, the accuracy is strongly de- 
graded at i + > 25.5. Table [J shows t he photometric 
redshifts measured by Illbert et al.l (|2009f h 

Most of the radio galaxies of our sample are in- 
cluded in the COSMOS Photometric Redshift Catalog 
Fall 2008, which provi d es the photometric redshifts mea- 
sured by Illbert et al] (|2009l ). We adopt again a 0'.'3 
search radius and find 29 objects. We report the re- 
sults in Table [3J One of these objects (236) turns 
out t o be a spectroscopic ally-confirmed Q SO at z = 
2.132jPrescott et al.ll2006h . The attempt of Illbert et all 
(|2009T l to estimate its photometric redshift failed because 
the templates used by these authors wer e not suitable to 
fit its AGN-dominated SED. Similarly. iMobasher etTal"] 
(|2007|) provided for this object a tentative photometric 
rcdshift, which is z=1.23±£^, which is clearly inconsis- 
tent with the spectroscopic value. 

6. RESUTS 

The SED modeling process has been performed for all 
the objects (except for the spectroscopically-confirmed 
QSO, object 236), by using the two template-fitting tech- 
niques (Fig ITS")) . As expected, 2SPD code has turned out 
to more reliably model the SEDs of our sources than Hy- 
perz. This is crucial in UV band. In fact, when significant 
UV emission appears in the SED, such as, in objects 31, 
258 and 285, Hyperz technique struggles to fit simulta- 
neously such a UV component and the remaining part 
of SED, which is well represented by an OSP. In this 
case, the best fit result from Hyperz is obtained using 
an old composite stellar population. This is due to the 
fact that the x 2 is lower if only the the old component of 
the SED is fitted, rather than the UV data points which 
are less well represented. Conversely, as expected, 2SPD 
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Fig. 6 . — Spectrum of ob j ect 25 (red solid line) from the Magellan 
survey (Trump ct al. 2007) which provides a spectroscopic redshift 
of 0.7917. We overlap the synthetic SED (green dashed line) de- 
rived from Hyperz at redshift z = 1.37, fitting the photometric 
points. The wavelengths on the top corresponds to the observed 
wavelengths, and at the bottom to rest frame. 

technique is more efficient to treat this case, since it just 
adds a small fraction of YSP to a dominant OSP to re- 
produce the UV emission, by departing from the concept 
of single star formation history. 

6.1. Photometric redshifts 

The first test on the accuracy of our photo-z derivation 
is a comparison with the spectroscopic redshifts, avail- 
able for 6 objects (namely 1, 16, 25, 31, 52, and 258). 
They are all compatible with each other within the er- 
rors with the only exception of object 25. For this ob- 
ject the photometric redshifts measured with Hyperz and 

2SPD (z = 1.37±g;?2 and z = L33 ±043' respectively) are 
significantly different from the redshift inferred from its 
Magellan spectrum (z = 0.7917). However, the Magellan 
spectrum docs not match with the COSMOS measure- 
ments, showing a large offset, of *~0.8 dcx. Apparently, 
the object observed by Magellan is not the radio galaxy 
25 and we do not consider its spectroscopic-z as reliable. 
Reassuringly, we checked that the spectra and photomet- 
ric data-points agree for the other 5 objects with available 
spectra. 

We then compare the photometric redshifts derived 
with the two SED fitting techniques. The overall range 
of z is in both cases from ^0.7 to ~ 3. The median 
photo-z are 1.21 and 1.30 for Hyperz and 2SPD respec- 
tively. Generally, the photometric redshifts derived with 
the two methods are consistent with each other within 
the errors (Fig. left panel). The normalized rcdshift 
differences (Az/(1 + z)) are smaller than 0.08 for all but 
5 objects that reach Az/z ~ 0.12 — 0.15., and a single 
strong outlier (object 37, Az/z ~ 0.28). For this galaxy, 
the Hyperz fit to the SED is particularly weak as it is not 
able to fit its photometric points at both FUV and MIR 
wavelengths. 

Furthermore, we compare the photometric redshifts 
measured wi th our method and with the template- fitting 
technique of Illbert etail (|2009D . Objects 2, 22, 28, 32, 
and 36 are not involved in this comparison since they are 
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TABLE 3 
Hyperz SED FITTING 
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Note. — Results from the analysis of the SEDs of the sample with Hyperz. Column description: (1) ID number of the object; (2) 
photometric redshift as result of Hyperz; all errors are quoted at at 99% of c onfidence level; (3) photom etric redshifts obtained with 2S PD; 
(4) source library of the best template which fits the observed SED: BC from lBruzual k. Ch ariot (20(H) and MA from lMar aston (2005D; (5) 
star formation history of the best template: ssp corresponds to a single stellar population and t = 'N' corresponds to the /^-models with 
exponentially decaying star formation rate with star formation timescale r = N Gyr; (6) age of the best template in Gyr; (7) Ay associated 
with the best template; (8) X^id °f ' ne best template; (9) stellar mass and its error of the galaxy in Log Mq as result of Hyperzmass. 



not included in the lllbert et all (|2009f) sample. Fig. [7] 
(right panel) shows the comparison between the photo - 
z resulting from 2SPD and from lllbert et all (|2009D . 
Again, the two methods yield similar photometric red- 
shifts with Az/(l+z) < 0.11 for all but 2 objects, namely 
226 and 228 for which we obtain Az/(1 + z) = 0.18 and 
0.25 respectively. In the case of object 226, similarly 
to the object 37 discussed above, the SED obtained with 
2SPD reveals a strong excess at both UV and MIR wave- 
lengths. Object 228 is instead a case of misidentification 
in the COSMOS multiband catalog. This result stresses 
the importance of our detailed work of identification to 
infer the genuine properties of the galaxies. 

Therefore, unless a spectroscopic redshift is available, 
we generally prefer to use the photometric redshift ob- 
tained with 2SPD because our code also includes the dust 
component which turns out to be important in the SED 
modeling. Table [5] shows the redshifts used throughout 



all the paper. 

6.2. Radio power distribution 

IChiaberge et al.l ([2009D selected the sample in the 
COSMOS field looking for radio sources with radio lu- 
minosities below the FR I/FR II break. Since we now 
have improved estimates of the photometric redshifts of 
the objects, we can derive a more accurate value for their 
radio power. We show in Fig. [8] the histograms of K- 
corrected (using a radio spectral index a = 0.8) radio 
powers at 1.4 GHz measured with NVSS and FIRST. 
The resulting 1.4 GHz luminosities are in the range 
1Q31-5 _ 1Q3 3-3 ergs -i Hz" 1 . The lo cal FR I/FR II break 
luminosity (jFanaroff fc Rilevlll974D . converted from 178 
MHz to 1.4 GHz adopting a = 0.8, is Li. 4GH z ~ 10 32 6 
erg s _1 Hz -1 . Therefore, the radio distribution of 
our sample actually straddles the FR I/FR II break. 
Nonetheless we remind the reader that radio sources with 
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TABLE 4 
2SPD SED FITTING 



ID 


rcdshift 


YSP 


OSP 


log JVL» 


Dust 


1R excess 


UV 




z phot 


Age 


A„ 
J\V 


Sysp 


IjOg 1V1* 


Age 






1 dust 


T , 


T r „ 


«8-24 


T 


1 


n OQ+U.U4 
U -° S -0.05 


u.uz 


u.ou 


±■0/0 


n zlq.% 
u.^±y /o 




n ir. 
U. i-O 


i n no + 0.04 


140 


2 4 








2 


i on + 0.10 

L - M -Q.Q9 


0.05 


0.90 


22.1% 




2.0 


0.98 


i i nn+0.04 

n.ooi ao4 


156 


7.4 


<43.65 




42.53 


3 


o on^O . 32 

2.20_ 044 


0.03 


0.12 


6 5% 


0.69% 


1.0 


0.00 


in Kn+0.08 
10.59_ . 10 


201-347 


250-59 


45.27 


-0.01 


43.21 


4 


i oy+0.10 

1,J '-0.06 


0.009 


1.00 


fi 9% 
u.y /o 


0.24% 


3.0 


0.57 


ii ic+0.04 
1L16 -0.03 


100 


5.2 


<43.55 




42.71 


5 


n n1 +0.22 
2 ' 01 -0.35 


0.02 


1.80 


53.1% 


5 q% 


3.0 


0.93 


11 ACI+0.04 

±±.4y_ 03 


121 


144.0 


44.76 


>0.49 




11 


1 cv+0.14 


0.007 


1.87 


1 7 9% 


0.62% 


3.0 


0.30 


m nc+0.10 


98 


15.4 


<43.80 






13 


i -m+0.08 


0.03 


1.17 


16.2% 


3 9% 


1.0 


0.58 


in -79+0.04 
10 - 72 -0.03 




7.0-6.2 


44.23 


-0.41 


42.65 


16 


n Q7+ 012 

UM ' -0.07 


0.006 


1.16 


14.8% 


0.23% 


2.0 


0.52 


m -7/1+0.06 
10 ' 74 -0.06 


158 


4.8 


43.60 


>-0.12 


iZAz" 1 


18 


n qo+0.14 


0.004 


2.55 


50.3% 


13.5% 


0.2 


1.35 


in n9+0 08 


116 


13.9 


<43.96 






20 


n oc+0.02 
U ' SB -0.02 


0.002 


1.26 


3 7% 
0. 1 /o 


0.07% 


3.0 


0.42 


11 n-j+0.02 
1L03 Io.03 


173 


1.4 


<43.02 




jn o f in 


22 


i Q n+0.05 


0.04 


1.18 


10.1% 


0.35% 


2.0 


1.50 


ii ic+0.02 
n ' 16 -0.03 


150 


11.8 


43.87 


>0.66 




25 


1 OO+0.11 


0.002 


1.17 


4.1% 


0.28% 


0.4 


1.05 


m 7^+0.04 


138 


8.8 


43.75 


>0.26 


in g TW, 
4Z.O / 


26 


1 no+ 012 

1 - ua -0.07 


0.006 


1.36 


7 4% 


n 1 7% 


1.0 


0.98 


ii 19+0.04 


128 


3.9 


<43.45 




42. 5U 


28 


9 Q n +0.20 
z.yu_ 26 


0.001 


1.97 


24.7% 


3.9% 


0.9 


0.00 


11 qc+0.04 
11.38_ 04 


146 


47.4 


44.28 


>0.22 




29 


i oq+0.23 
1 -' 3z -0.24 


0.004 


0.38 


11.9% 


0.18% 


0.4 


0.94 


in m+0 05 


83 


16.4 


<43.68 




42.83 


30 


L06 -0.07 


0.001 


1.77 


3 n% 

o.u /o 


0.13% 


2.0 


0.73 


11 nQ+0° 5 
1L03 -0.05 


130 


3.5 


<43.47 






31 


n oq+0.03 
U - 8S -0.05 


0.007 


0.23 


7 5% 

( .0/0 


0.06% 


2.0 


0.23 


1n -,..+0.03 
10 ' 75 -0.03 


164 


2.6 


<43.35 




42.86 


32 


9 71 +0.38 
z - ' -0.34 


0.002 


0.98 


10.6% 


0.26% 


2.0 


0.39 


in no+0.04 

10.98r (J . O4 




51-58 


44.96 


-0.98 


43.05 


34 


1 cc+0.41 
1 - oa _0.19 


0.005 


0.55 


11.0% 


0.05% 


3.0 


0.51 


i n nn+0 07 
iu.yy_ 07 


140 


19.7 


<44.04 




42.82 


36 


i nv+0'10 
1 - 07 -0.04 


0.008 


2.30 


42.4% 


2.0% 


3.0 


1.07 


in oq+0.02 
10.83_ 03 


186 


4.4 


43.46 


>0.20 




o t 


1 -J°_0.42 


0.001 


0.00 


1 n% 


0.09% 


0.3 


0.09 


in fil+0.11 


213-614 


84-22 


45.03 


-0.25 


44 10 


38 


i n n +0.17 

i.ou_ 2g 


0.005 


0.95 


27.0% 


0.58% 


0.9 


0.86 


m ck+0.07 

10.65I 007 


260 


6.2 


43.77 


0.45 


43.03 


39 


i 1n +0.05 
L10 -0.05 


0.007 


1.86 


15.0% 


0.83% 


1.0 


0.97 


in ac+0.03 
10.88_o.o3 


80 


5.3 


<43.20 






52 


n 7/ ,+0.18 
°" 4 -0.19 


0.004 


0.66 


13.6% 


0.11% 


2.0 


0.34 


in -7S+0.10 
10 ' 78 -0.10 


189 


2.6 


43.33 


0.61 


42.96 


70 


£i ofi^O . 53 
-0.20 


0.009 


0.48 


12.8% 


0.87% 


0.4 


0.74 


10 65+ ' 07 
lU.DO_o 03 


180 


27.5 


44.07 


>-0.50 


43.54 


202 


1 O1+0.09 

i " 31 -0.12 


0.005 


0.14 


1.2% 


0.004% 


3.0 


0.14 


10.86t°$ 


124 


7.4 


<43.66 




42.21 


219 


i n3 +0 02 

1 - Uci -0.04 


0.02 


0.83 


6.2% 


0.69% 


0.3 


1.47 


10 67+ ' 03 


112 


2.8 


<43.32 




42.69 m 


224 


1 10+ 010 
i - ±u -0.04 


0.004 


0.63 


3.8% 


0.04% 


1.0 


0.76 


10 7 i+0.03 
±u -' 1 -0.04 


155 


4.5 


<43.53 




42.53 m 


226 


9 ocr+0.63 
i,OJ -0.31 


0.001 


0.23 


2.8% 


0.12% 


0.9 


0.00 


10 5Q+ ' 07 
±u.oy_ og 


212-467 


88-36 


45.01 


-0.46 


43.34 


228 


1 O1+0.05 


0.002 


2.29 


37.1% 


0.72% 


2.0 


2.21 




133 


6.9 


<43.64 






234 


1 10+ ' 14 


0.006 


2.30 


46.3% 


2.2% 


2.0 


1.05 


io.83+8:8 


104 


5.8 


<43.52 






258 


96+ 19 
u ' yD -0.13 


0.001 


0.08 


1.6% 


0.05% 


0.08 


1.72 


10 64+ 08 


134 


3.0 


<43.41 




43.27 


285 


1 10+ 013 
i - ±u -0.08 


0.001 


0.13 


2.1% 


0.05% 


0.04 


1.56 


in 4 q+0-04 
±u.4o_ 0(H 


135-448 


1.9-1.8 


43.73 


-0.22 


43.08 



Note. — Results from the analysis of the SEDs with 2SPD. Column description: (1) ID number of the object; (2) photometric redshift 
measured with 2SPD; (3)-(4)-(5)-(6) age in Gyr, Ay, flux fraction and mass fraction of the young stellar population (YSP) at 4800 Arest 
frame; (7)-(8) age in Gyr and Ay of the old stellar population (OSP); (9) the total stellar mass of the galaxy in Mq; (10)-(11) the effective 
temperature (in K) of the one or two dust components and their luminosities, l^dust ( m units of 10 9 L©); (12)-(13) the infrared excess 
luminosity (in erg s _1 ) defined in the text (Section 16.41 1 and the spectral index measured on the infrared excess at 8 and 24/^m; (14) UV 
luminosity at 2000 A in the rest frame in erg s — 1 . The marginal UV excesses are marked with a m . 



a clear FR II morpholo gies were previously excluded by 
IChiaberge et~aT1 (|2009l ). 

As already discussed by IChiaberge et al.l (|2009D . the 
high-frequency radio observations available for this sam- 
ple (their rest frame frequency are in the ^3-5 GHz 
range) are not ideal to observe and detect their total ra- 
dio structures because they miss diffuse/steep-spectrum 
radio emission. Lower frequency radio observations are 
needed for a more accurate measurement of their intrin- 
sic total radio power distribution. 

Since the radio power distribution is rather broad (cov- 
ering two orders of magnitudes) , we separate the sample 
into two groups, including sources of high and low power 
(HP and LP, respectively) using the local FR I/FR II 
break as divide (see Table [SJ . This operative definition 
enables us to explore the role of radio power in defining 
the overall properties of these galaxies. Two third of the 
sample shows radio luminosities below the break. The LP 



and HP classes are, not surprisingly, roughly separated 
also in redshifts (Fig. since our sample is limited in 
radio flux, and indeed the HP sources generally lie at z 
> 1-5- 

6.3. Host galaxy properties inferred from SED modeling 

We now focus on the properties of the host galaxies 
inferred from SED modeling and in particular on their 
stellar populations. 

One of the key and more robust result is the mea- 
surement of the stellar content, M*. The inferred mass 
range is ~10 10 — 10 12 M Q . The masses measured with 
Hyperzmass and with 2SPD are compared in Fig. 1101 
The two techniques return generally consistent values of 
the stellar masses for most sources (the median values 
are in both cases ~ 7 x 10 10 M Q ), but with a few ev- 
ident outliers. The stellar masses derived from singlc- 
and two- stellar components are know to differ signif- 
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Fig. 8. — Distribution of the K-corrected total radio luminosity 
(in erg s _1 Hz -1 ) at 1.4 GHz from the FIRST data (upper panel) 
and NVSS data (lower panel). The dashed lines corresponds to the 
local FR I/FR II break used to separate the sample in high power 
(HP, shaded histograms) and low power (LP) sources. 



icantly (e.g. IPapovich et al.l [20061 ) . Furthermore, the 
inclusion of dust components by 2SPD has two effects 
by changing i) the stellar mass since thermal emission 
largely contributes to observed infrared emission, and ii) 
the photometric redshift. The largest outlier is object 37 
whose SED is poorly fitted by Hyperz with a very high 
mass 7.0 x 10 11 M Q , while a strong dust component and 
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Fig. 9. — Rcdshifts versus the K-corrected FIRST radio luminos- 
ity at 1.4 GHz (in erg s~ x Hz" 1 ). The empty points are the LPs 
and the full points are the HPs. 

a mass 10 times lower are required by 2SPD. At the op- 
posite end of the mass distribution, we find two objects 
(36 and 226) in which again 2SPD finds a strong dust 
component, but where the Hyperz is a factor 3 to 7 lower 
than that obtained with 2SPD. Apparently, the presence 
of dust components has a stronger impact on the stel- 
lar mass estimate than the inclusion of a young stellar 
population. 

Considering the radio power of the sources (see Fig.fTTl 



12 



Baldi R. D. et al. 



12 pr 

10 E- 

8 - 

2? 6 E- 

4 E- 



n i i i r 




10 10.5 11 11.5 12 
Log M, (2SPD) 



12 r 

10 E 
8 E 
6 E 
4 E 
2 E 
E 



T 



t I i i r 



12 



11.5 
1 1 

w 11 



10.5 




10 



* 




M — i 



10 



10.5 11 11.5 
Log M» (Hyperz) 



1 2 



10 10.5 11 11.5 12 
M, (Hyperz) 



Fig. 10. — Left panel: distribution of stellar masses (in Mq) of our sample obtained with 2SPD (upper panel) and with Hyperz (lower 
panel). Right panel: comparison between the stellar masses (in Mq) inferred from SED modeling using Hyperz and 2SPD techniques. The 
black histogram (points) represents the HPs, while the empty one the LPs. 



left panel), HPs show slightly larger stellar mass con- 
tent than LPs. The median values (obtained from 2SPD 
code) for LPs is 5.8 x 10 10 M Q , while for HPs is 8.7 x 
10 10 M Q . However, the two distributions are very broad 
and show an almost complete overlap. A very similar 
result is found looking for differences in mass between 
objects at different redshifts (see Fig. [TTJ right panel). 

Due to the degeneration inherent to the modeling of 
SEDs the remaining results of the stellar populations are 
relatively unconstrained. Overall, if we concentrate in 
the whole sample, the SED fitting shows that the hosts 
are mainly dominated by a massive OSP with an age of 
~ 1—3 x 10 9 years. Although less reliable, the YSPs 
required to fit the UV excesses show ages of 1 — 30 Myr, 
with a contribution to the total mass of the galaxy of 
< 1% for most sources. The flux contribution at 4800 A 
(rest frame) of the YSP is <30% for most of the objects. 

6.4. Dust emission 



As discussed in Section 15.21 the parameters related 
to the dust emission must be taken as approximate 
rather than a real measurements of the dust component. 
Nonetheless, it is still significant that dust emission is re- 
quired to adequately model the SEDs of 15 objects due 
to the detection of emission at 24 fim. In addition sig- 
nificant excesses above the stellar emission are observed 
also at shorter infrared wavelengths in 8 of these galaxies. 
In order to explore the dust properties we estimated the 
residuals between the best fitting model including only 
the stellar component and the data-points, looking for 
an excess at long wavelengths due to dust. We then inte- 
grated the residuals to obtain the total dust luminosity 
in the range covered by the Spitzer data, i.e. ~ 3 - 26 
fim. The integration is performed by assuming that the 
spectrum is represented by a multiple step function (see 



Fig. EH. 

The resulting dust luminosities (see Table 0J, esti- 
mated as IR excess luminosities, are in the range Ldust ~ 
10 43 - 5 - 10 45 - 5 erg s -1 . A trend links radio and dust 
emission (see Fig. [T3l left panel) with most (7/9) HP 
radio-sources showing a significant dust emission with 
luminosities larger than L^ust ~ 10 44 erg s _1 . The LPs 
are instead (with only 2 exceptions) of lower Ldust and 
in many cases (16 galaxies) only upper limits can be de- 
rived. 

We check the statistical significance of such a trend, 
by using a censo red statistical analysis (ASURV, 
lLavallev et al.lfl992T ) which takes into account of the pres- 
ence ofjap_perjimits. Using the generalized Kendall's r 
test (Kendall 1983), the probability that a fortuitous cor- 
relation appears is 0.0003. However, the common depen- 
dence of the two luminosities on redshift might play an 
important role. Therefore, we perform a partial corre- 
lation analysis (jAkritas Sz Siebertlll996t) to examine the 
linear relation between the luminosities excluding the de- 
pendence on redshift. Opcratively, we use the partial 
Kendall's r test, whose null hypothesis is the absence 
of the correlation excluding the redshift variable. The 
partial Kendall's coefficient r is 0.20 and the standard 
deviation is 0.081. The null hypothesis is then rejected 
at the level of 0.05. 

We also measured the spectral index of the residuals 
between 8 /zm and 24 /im, Q!(8_24)- Taking into account 
only significant excesses at 8 /im (residuals larger than 
3 a), this value can be estimated in 8 cases, with values 
spanning between a(g_24) ~ 1 and -1. For the remaining 
objects with only a 24 /im detection, the limit to the 8 
/im flux translates into a lower limit of ct( S _24\ > 1. This 
parameter can be crudely related to the overall dust tem- 
perature. By assuming a single black-body dust compo- 



Study of SEDs of radio galaxies at z~d-3 



13 



Q 

a* 

m 

C\2 



11.5 



11 



60 10.5 



o 



10 



n i i r 



o 

o o 



31.5 



o o 
P o 



_l I I L 



C\2 



11.5 



11 



60 10.5 



o 



10 



oo 
o 





om 


o 







o 


o 


o 

CO 


o 


o°° 


o 




• 




o 



o o o 



32 32.5 
Lop L 



33 



33.5 



0.5 



1.4 GHz. FIRST 



1.5 2 2.5 
redshift z 



Fig. 11. — Stellar masses (in Mq) measured with 2SPD in relation with the K-corrected FIRST radio powers (in erg s 1 Hz 1 ) (left 
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Fig. 12. — Residuals obtained by subtracting the stellar compo- 
nent from the SED in the 5-bands Spitzer data which we ascribe to 
the dust emission. We show the results for object 3 (black) and 5 
(red) as an example. The step function connecting the data-points 
is used to obtain the IR excess flux, i.e., the total dust flux. 



ncnt, the values of Q!(g-24) translate into a temperature 



range of 500-850 K and 300-550 K for a 



(8-24) 



-1 and 



1, respectively. The derived temperature depends on red- 
shift, with the lower (upper) values of T being derived 
for z = 0.75 (z = 2). 

Fig. [T3] (middle panel) shows a broad relation between 
the IR spectral index and the IR excess luminosities. 
This would mean that the decrease of «(s-24) is ascrib- 
able to the increase of the high-temperature dust compo- 
nent. Since a large radio power also implies a large dust 
luminosity, this may suggest a possible AGN nature of 
the dust heating source for the objects with largest IR 



6.5. UV excess 



excess. 



As noted in Sect. 15.11 the Hyper z code often does not 
reproduce satisfactorily the bluest part of the SED and 
we chose to model such part of spectrum with a young 
stellar population, included by our 2SPD. The inclusion 
of this component does not alter significantly the photo- 
metric redshift or the galaxies mass, but it is clearly of 
great importance to understand the nature of these radio 
galaxies. Nonetheless, inspection of the SEDs obtained 
with 2SPD indicates that the UV excesses (above the 
contribution of the old stellar component) are usually 
very poorly constrained. Furthermore, the very stellar 
origin is not granted and the UV excess might hide an 
AGN contribution. 

It is necessary to introduce a model-independent crite- 
rion to assess which sources actually show an UV excess 
and to estimate its luminosity. We visually inspected all 
SEDs, searching for sources with a substantial flattening 
in the SED at short wavelengths or with a change of the 
slope between the OSP and the emission in the UV band. 
Fifteen sources show a clear UV excess (namely object 
2, 3, 4, 29, 31, 32, 34, 37, 38, 52, 70, 202, 226, 258, and 
285) and additional 7 objects show a marginal UV excess 
(namely object 13, 16, 20, 25, 26, 219, and 224). In the 
remaining galaxies the SEDs in UV bands drop sharply 
and are well reproduced by the emission from OSPs. In 
order to quantify the UV contribution, we measure the 
flux at 2000 A in the rest frame, Lyy, from the best 
fitting model, for those objects showing the UV excess. 
The UV luminosities range is 10 42 — 10 44 erg s _1 . HPs 
show larger UV luminosities than LPs by a factor 2.5, 
even though the strongest UV excess is seen in the LP 
radio galaxy 37. If we concentrate the sources which do 
not show an UV excess, the 2000 A luminosity is < 10 43 



14 



Baldi R. D. et al. 



45 - 



45 - 



o 



£44 



43 - 



40 




o 
-J 



£44 



43 - 



1 


III 


II 


III 

• 


1 1 1 1 1 




— 


• 


• 


o 






_ 












- 









•> 


- 


- 




•» 






- 










G» 





































o> 












O 




1 


III 


III 


III 


1 1 1 1 1 1 





41 42 

Lj 4 GHZi FIRST 



43 



1.5 



-0.5 



0.5 



a 



IK 
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spectral index from 8 to 24 (im estimated from the IR excess in the SED above the stellar emission. Empty points are the LPs and full 
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erg s _1 (mostly, ~ 10 42 erg s^ 1 ). 

In order to explore the nature of the UV emission, it is 
useful to compare the UV luminosity with the multiband 
properties obtained in the previous sections. No relations 
appear to link the UV excess to the redshift and the mass 
of the galaxy. Conversely, the UV luminosity appears to 
be linked to the radio power and IR excess luminosity 
(Fig. ED. 

For the radio and UV luminosities the linear correla- 
tion coefficient is r = 0.404 and the probability of the 
presence of a fortuitous relation is P= 0.062. The partial 
Kendall's coefficient r is 0.17 and the standard devia- 
tion is 0.15: the null hypothesis, i.e. the absence of the 
correlation, can not be rejected. 

For the IR-UV relation, the generalized Kendall's r 
test returns a probability of no correlation of P= 0.0072. 
We obtain a partial Kendall's coefficient r = 0.41 and a 
standard deviation of 0.069: the null hypothesis, i.e. the 
absence of the correlation, is rejected at the level of 0.05. 

7. DISCUSSION AND CONCLUSIONS 

We used the multiwavelength data provided by the 
COSMOS survey to analyze the first sizeable sample of 
34 low-p ower radio galaxies at h igh redshifts (z > 1), se- 
lected bv lChiaberge et al.l (|2009l ). We performed a care- 
ful visual inspection of all their multiwavelength counter- 
parts to identify their genuine emission and to infer their 
SEDs. Those would have been compromised in case of a 
blind use of the COSMOS catalog. Taking advantage of 
the FUV-MIR spectral coverage, we modeled their SEDs 
using two different template fitting techniques: i) the 
Hyperz code that uses a single stellar population and ii) 
our own code 2SPD that includes also dust component (s) 
and a young stellar population. We analyzed the proper- 
ties of the SEDs of the sample. We here summarize the 
main results and briefly discuss them: 



(i) The photometric redshifts of these radio sources 
range from ~0.7 to 3. The photo-z measured with 
the two techniques are generally con sistent with 
each other, with those measured by lllbert et"aLl 
(|2009l ). and with the available spectroscopic red- 
shifts. In addition, we measured the photometric 
red shifts for five objec ts which were not included in 
the lllbert et~aTI (|2009f ) sample. 

(ii) The new and accurate measurements of the pho- 
tometric redshifts enable us to infer the radio 
power distribution of these objects. The result- 
ing K-corrected 1.4 GHz luminosities are in the 
range 10 315 - 10 33 3 erg s _1 Hz -1 , straddling the 



FR I/FR II break luminosity (Li. 4G h z ~ 10 32 6 erg 
s _1 Hz -1 ) as defined for low redshifts galaxies. 



(iii) The resulting stellar masses are mostly confined in 
the range ^io 10 - 5-11 - 5 M Q , with a median value of 
7 x 10 rt Mq. 



(iv) The SED of most objects of the sample is consistent 
with the presence of a dominant contribution from 
an old stellar population with an age ~ 1 — 3x 10 9 
years. However, significant excesses are often ob- 
served at the shortest and/or longest wavelengths. 

(v) A dust component is needed to account for the 
24 fim emission in 15 galaxies and significant ex- 
cesses above the stellar emission are observed also 
at shorter infrared wavelengths in 8 of these galax- 
ies. Estimates of the dust luminosity yield values 
in the range L dnst ~ 10 43 5 - 10 45 5 erg s -1 . The 
overall dust temperature, estimated for the 8 radio- 
galaxies with a substantial dust excess at A < 8fim, 
is in the range 300-850 K. 
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Fig. 14. — UV luminosity (ergs 1 ) measured at 2000 A versus: (left panel) K-corrected FIRST radio luminosity (ergs 1 ) , and (right 
panel) infrared excess luminosity ( ergs" 1 ) . Empty points are the LPs and full points are the HPs. 



(vi) 



(vii) 



Inspection of the SED obtained with 2SPD indi- 
cates that the UV excesses (above the contribution 
of the old stellar component) are often present (sig- 
nificantly in 15 sources and marginally in 7 sources), 
but they are usually weakly constrained. The UV 
luminosities measured at 2000 A (rest frame) is in 
the range 10 42 — 10 44 erg s _1 . 

Although the censored analysis does not provide 
significantly high statistical parameters, we can ten- 
tatively confirm the presence of positive links be- 
tween the dust emission with both the radio and 
UV luminosities. For these relations, the possibility 
of a common luminosity dependence on distances is 
rejected at the level of 0.05. 



The selection performed by IChiaberge et al.l (|2009f ) 
with the aim of searching for FR I candidates at z ~ 1 — 2 
turned out to be successful. In fact, our work confirms 
i ) their location in the range of redshifts aimed with 
the selection process, although extending slightly below 
1 and up to z ^ 3 and ii) the low radio luminosity of the 
sample, generally consistent with those of the local FR I, 
although 1/3 of the sample exceeds the local FR I/FR II 
luminosity break. The extension of their radio power to 
larger luminosities does not necessarily imply a FR II 
nature for those sources. In fact the sources with a clear 
FR II morphology were excluded from the analysis. In 
addition, in the local Universe, the radio distribution of 
FR I s is broad and overla ps with the FR II distribu- 
tion (jZirbel k Bauml IT99 5f) . Furthermore, this overlap 
increases at higher radio frequencies. 

Overall, the hosts of these high-z low-luminosity radio- 
sources are similar to those of the local FR I whic h usu- 
ally l i ve in red massive early-type galax i es (e. g .. iZirbel 
1996 1 iBest et all [2005] IBaldi k CapettH 120081 iSmolcic 
2009| IBaldi k CapettfelOD . 



However, the SED modeling reveals that additional 
components to the old stellar population have to be in- 
cluded to account for the emission at the SED extremes, 
i.e. in either the UV or in the MIR band, in most of 
the sources of the sample. This behavior is not seen in 
the low-z FR I that are generally faint in UV (both from 
the point of view of star formation and nuclear emis- 
sion, lChiaberieiFin[20ni |BaId[Sciiitl[200l). Sim- 
ilarly, their MIR luminosities exceed by a very large fac- 
tor (between 30 and 3000 for the MIR detected sources) 
the typical low-z FR I luminosities (~ 10 42 ergs -1 , 
iHardcastle et all 12009( 1. Conversely, the UV and MIR 
properties are somewhat similar to those of local FR lis , 
which show bluer color (e.g., IBaldi k Capettil 120081: 
Smolcidl2009f) and large dust amount (e.g. Jde Koff et all 
2000t iDicken et al.ll2010l) than FR Is. 

The origin of the MIR and UV emission can not 
be firmly established based on the available data. 
The estimate of the dust temperature (possible for 
only 8 objects) is in the range expected from dust 
heati ng from a quasar-l i ke n u cleus (T ^ 30 K, 
e.g. iSiebenmorgen et al.l 120041: I Ogle et al.M2006D and 
far larger than the dust associated with star forma- 
tion (e.g.. iHwang et al.ll2010t iSreenilavam k Fichl I20T1 



iBoquien et al.ll2011b iPatel et al.ll2011h . The high end of 
the MIR luminosity reaches values similar to those found 
in high power radio-galaxies and QSOs. However, dust 
emission is seen only in less than half the sources of the 
sample. 

Similarly, a UV component in excess to the old stel- 
lar population is firmly detected in the same number of 
sources (but there is not a one-to-one correspondence 
between UV and MIR emission). The observed UV lu- 
minosities are much larger that the faint non-thermal 
UV-nuc lei seen in FR I and, instea d, similar to those 
FR lis (jChiaberge et al.lll999l I2002D . dominated by the 
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TABLE 5 

'Final' redshifts and radio properties of the sample 



ID 


z 


-Ljvvss 


Lfir.st 


class 


1 


0.88 s 


<31.78 


31.93 


LP 


2 


1.33 


32.02 


32.4 


LP 


3 


2.20 


33.1 


33.19 


HP 


4 


1.37 


32.77 


32.86 


HP 


5 


2.01 


32.47 


32.89 


HP 


11 


1.57 


<32.18 


32.52 


LP 


13 


1.19 


32.02 


32.22 


LP 


16 


0.97 s 


32.38 


32.27 


LP 


18 


0.92 


32.22 


32.28 


LP 


20 


0.88 


<31.66 


31.93 


LP 


22 


1.30 


<32.38 


32.34 


LP 


25 


1.33 


32.29 


32.39 


LP 


26 


1.09 


32.02 


32.25 


LP 


28 


2.90 


32.99 


33.13 


HP 


29 


1.32 


32.28 


32.31 


LP 


30 


1.06 


31.83 


32.11 


LP 


31 


0.91 s 


32.14 


32.18 


LP 


32 


2.71 


32.8 


33.17 


HP 


34 


1.55 


32.84 


32.77 


HP 


36 


1.07 


32.23 


32.25 


LP 


37 


1.38 


32.26 


32.33 


LP 


38 


1.30 


32.94 


33.00 


HP 


39 


1.10 


<31.90 


32.16 


LP 


52 


0.74 s 


<31.54 


31.73 


LP 


70 


2.32 


33.12 


33.18 


HP 


202 


1.31 


31.98 


32.46 


LP 


219 


1.03 


<31.96 


32.09 


LP 


224 


1.10 


32.28 


32.27 


LP 


226 


2.35 


<32.61 


32.94 


HP 


228 


1.31 


32.25 


32.51 


LP 


234 


1.10 


32.41 


32.48 


LP 


236 


2.13 


33.29 


33.29 


HP 


258 


0.90 s 


31.90 


32.12 


LP 


285 


1.10 


32.23 


32.31 


LP 



Note. - Column description: (1) ID number of the ob- 
ject; (2) redshift of the object used throughout the work, 
with the spectroscopic redshifts marked with a s ; (3)-(4) 
K-corrected radio luminoisty at 1.4 GHz from NVSS (from 
http://www.cv.nrao.edu/nvss/NVSSlist.shtml) and FIRST (taken 
from IChiabcrgc et al.l I2009T) in erg s 1 Hz 1 ; (5) classification 
based on the radio power: low or high power (LP or HP) radio 
sources. 

accretion disk emission. 

Summing up, we find that the sources of the sample 
display a wide variety of properties, despite the relatively 
narrow range in radio luminosities. The 8 objects with 
the strongest MIR excess, and with high dust tempera- 
ture indicative of a quasar-like nature, also show a sig- 
nificant UV excess. To this group we must obviously add 
the object 236, the spectroscopically confirmed QSO. We 
must note, however, that the SED of these 8 sources is 
much redder than that of object 236 and that their HST 
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images do not provide any evidence fo r the presence of 
bright unresolved nuclei, as described bv lChiaberge et al] 
(|2009f) based on visual inspection of the data. 

At the opposite end of the radio luminosity distribu- 
tion, there are 7 galaxies not showing any UV excess nor 
any 24 /xm emission (and this number raises to 11 by in- 
cluding also the objects with only marginal UV excesses) . 

The remaining galaxies, amounting to about half of 
the sample, show an excess only in either the UV or MIR 
band. For those objects which do not show an UV excess 
or arc not detected at 24 /zm, the limits to the UV/MIR 
ratio are broadly consistent with the values obtained for 
the remaining galaxies of the sample. In particular we 
do not find objects with a high MIR luminosity without 
an UV excess that might be expected in the case of an 
obscured QSO. Unfortunately, with the available infor- 
mation we cannot distinguish between an origin related 
to star formation or to an active nucleus. In addition, 
relatively large amounts of dust, suggested by the dust 
luminosities, indicate that possible obscuration may pre- 
vent us from detecting UV emission. 

A further detailed analysis of the nuclear properties is 
needed to understand which type of AGN are associated 
with these high- z radio galaxies. T his will include their 
X-ray emission ()Tundo et al.l 1201 1[ ) and the radio core 
flux, available from the COSMOS/VLA data, but which 
we defer to a future study. This might also provide new 
insights on the con t roversial FR I-QSO association (e.g., 
Falcke et al.l [19951: iBaum et al.l [l99ol iCao k Rawlings 
2004t iBlundell k Rawlingsl 120011 see also iBlundelll 12003 
for review on this subject). 

Summarizing, this work validates the first sizeable 
sample of low-luminosity radio galaxies at high redshifts, 
0-7 z < 3. This opens the possibility to perform a 
detailed comparison of the host and nuclear properties 
of these sources with those of i) the local low-luminosity 
radio-galaxies, ii) the powerful radio-galaxies in the same 
redshift range, and iii) the population of non active 
galaxies at z ~ 1 — 3. These issues will be addressed 
in a forthcoming paper. 
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Fig. 15. — Left panels: SEDs of the sample (green line) fitting the photometric points, as result of Hyperz. We plot the photometric 
point at 24 fim (Spitzcr/MIPS), if available, but is not included in the modeling. Right panels: SEDs of the objects fitted with 2SPD. The 
total model is the green line, the YSP is the blue line, the OSP is the red line, and the dust component(s) is the light blue line. For both 
the panels, the wavelengths on top of the plots correspond to observed wavelengths, while those on bottom are at rest frame. For object 
236 (QSO) we only show the photometric points 



TABLE 6 

COSMOS MULTIWAVELENGTH COUNTERPARTS OF THE SAMPLE 





ID 


1 

u* 


Bj 


3 + 


Vj 


r+ 


i* 


i + 


F814W 


z + 




J 




K S 




K 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




(11) 




(12) 




(13) 




1 


26.56±0.15 


25.89±0.13 


25.61±0.10 


24.42±0.05 


23.55±0.03 


22.24±0.05 


22.29±0.11* 


21.87±0.06 


21.32±0.01 


20 


46±0 


10* 


19.59±0.09* 


19 


57±0 


03 


2 


26.79±0.20* 


26.43±0.11* 


26.29±0.11* 


25.93±0.10* 


25.26±0.08* 


24.35±0.45* 


24.43±0.10* 


24.49±0.34* 


23 . 44 ± . 07 * 


22 


18±0 


15* 


20.85±0.02 


20 


90±0 


15 


3 


26.69±0.17 


25.58±0.10 


25.81±0.12 


25.60±0.11 


25.53±0.10 


25.31±0.78 


25.10±0.08 


24.96±0.79 


24.48±0.12 


23 


39±0 


25* 


21.57±0.03 


22 


84±0 


96 


4 


26.64±0.17 


25.93±0.13 


25.88±0.13 


25.46±0.10 


25.12±0.08 


24.27±0.36 


24.20±0.05 


23.82±0.29 


23.27±0.04 


21 


93±0 


10* 


20.53±0.01 


20 


74±0 


16 


5 


<25.77* 


<25.68* 


<25.53* 


<25.17* 


<25.23* 


24.12±0.27 


24.59±0.41* 


23.98±0.33 


23.41±0.25 


22 


24±0 


13* 


20.65±0.12 


20 


82±0 


15 


11 


<29.40* 


28.75±1.14 


28.00±0.76 


27.11±0.35 


26.27±0.17 


25.16±0.84 


24.96±0.07 


24.64±0.60 


24.09±0.08 


22 


28±0 


07* 


21.18±0.16 


21 


75±0 


49 


13 


26.14±0.15* 


25.75±0.13* 


25.74±0.13 


25.02±0.09* 


24.38±0.06* 


23.20±0.17 


23.29±0.04* 


22.90±0.13 


22.21±0.02 


21 


44±0 


09* 


20.46±0.10 


20 


54±0 


16 


16 


26.07±0.12 


25.71±0.12 


25.32±0.11 


24.84±0.08 


24.23±0.05 


23.10±0.14 


23.15±0.03 


22.85±0.10 


22.16±0.02 


21 


47±0 


08* 


20.40±0.08 


20 


50±0 


11 


18 


26.16±0.13* 


25.76±0.10* 


25.51±0.10 


24.61±0.08* 


23.99±0.06* 


22.91±0.11* 


22.80±0.04* 


22-51±0.10* 


22.04±0.03* 


21 


54±0 


07* 


20.71±0.04 


20 


92±0 


18 


20 


26.11±0.15 


25.74±0.15 


25.51±0.10 


24.57±0.06 


23.67±0.04 


22.35±0.09 


22.40±0.02 


22.19±0.08 


21.53±0.02 


20 


79±0 


06* 


19.91±0.07 


20 


10±0 


07 


22 


<27.38* 


<27.44* 


<26.97* 


<26.85* 


<26.00* 




<25.13* 




23.93±0.06 


22 


14±0 


17* 


20.79±0.10 


21 


24±0 


15 


25 


26.66±0.22 


25.36±0.09 


25.44±0.10 


24.97±0.07 


24.55±0.05 


23.64±0.32 


23.52±0.03 


23.32±0.19 


22.57±0.03 


21 


43±0 


09* 


20.57±0.05 


20 


59±0 


17 


26 


26.58±0.20 


25.44±0.10 


25.56±0.10 


24.92±0.07 


23.90±0.04 


22.68±0.10 


22.60±0.02 


22.32±0.09 


21.57±0.01 


20 


42±0 


04* 


19.42±0.03 


19 


48±0 


04 


28 


<27.39* 


<27.50* 


<27.20* 


<26.90* 


<26.80* 




25.55±0.30* 




24.27±0.20* 


22 


48±0 


20* 


20.99±0.08* 








29 


25.81±0.12 


25.51±0.10 


25.49±0.09 


25.31±0.09 


25.18±0.08 


25. 00 ±0.87 


24.78±0.07 


24.92±0.74 


24.30±0.11 


22 


86±0 


27* 


22.19±0.10 


22 


89±0 


80 


30 


27.85±0.63 


26.77±0.22 


26.82±0.23 


25.90±0.13 


24.69±0.06 


23.41±0.19 


23.41±0.03 


23.01±0.15 


22.32±0.02 


21 


25±0 


08* 


20.02±0.04 


20 


19±0 


09 


31 


24.75±0.08 


24.50±0.0e 


24.51±0.05 


23.94±0.04 


23.42±0.03 


22.33±0.08 


22.33±0.02 


22.07±0.07 


21.61±0.02 


20 


97±0 


08* 


20.12±0.05 


20 


32±0 


09 


32 


<27.54* 


27.85±0.48* 


<27.19* 


<27.00* 


26.44±0.30* 




26.07±0.25* 




25.15±0.28* 


24 


13±0 


50* 


22.43±0.17* 


23 


24±1 


44 


34 


26.09±0.40* 


26.19±0.20* 


26.00±0.19* 


25.58±0.12* 


25.44±0.10* 


25.09±0.66* 


24.88±0.07* 


24.95±0.25* 


24.12±0.11* 


23 


10±0 


51* 


21.10±0.29* 


20 


96±0 


14 


36 


<27.8* 


<26.80* 


<26.40* 


26.18±0.45* 


25.50±0.31* 


24.31±0.33* 


24.15±0.15* 


24.00±0.25* 


23.20±0.24* 


22 


04±0 


15* 


20.62±0.07 


20 


44±0 


10 


37 


22.70±0.02 


22.35±0.02 


22.63±0.02 


22.38±0.02 


22.18±0.02 


21.93±0.05 


21.93±0.01 




21.17±0.01 


20 


31±0 


05* 


20.14±0.11* 


20 


06±0 


08 


38 


25.50±0.10 


25.00±0.08 


24.96±0.07 


24.66±0.06 


24.20±0.04 


23.58±0.32 


23.49±0.03 


23.47±0.22 


22.69±0.03 


22 


08±0 


16* 


20.89±0.05 


21 


34±0 


22 


39 


<27.07* 


26.05±0.15* 


<26.00* 


25.28±0.09 


24.42±0.05 


23.13±0.18 


23.13±0.03 


22.86±0.15 


22.18±0.02 


21 


06±0 


08* 


19.99±0.04 


20 


27±0 


09 


52 


24.24±0.04 


23.64±0.04 


23.76±0.04 


23.13±0.03 


22.62±0.02 


21.62±0.05 


21.67±0.01 


21.52±0.05 


21.18±0.01 


20 


58±0 


05* 


19.79±0.04 


19 


78±0 


06 


70 


26.16±0.14* 


25.14±0.09* 


25.27±0.14* 


24.79±0.08* 


24.51±0.06 


24.34±0.58 


24.35±0.06 


24.13±0.43 


24.01±0.09 


2 3 


42±0 


50* 


21.76±0.09 


21 


40±0 


27 


202 


27.03±0.28 


27.11±0.29 


27.40±0.42 


26.92±0.29 


25.86±0.12 


24.44±0.53 


24.41±0.06 


24.29±0.45 


23.48±0.05 


22 


19±0 


17* 


21.07±0.05 


21 


62±0 


30 


219 


25.68±0.10 


25.09±0.08 


25.23±0.08 


24.49±0.05 


23.81±0.04 


22.74±0.09 


22.71±0.02 


22.37±0.09 


21.70±0.02 


21 


02±0 


06* 


19.93±0.04 


20 


15±0 


37 


224 


26.40±0.18 


25.75±0.13 


25.52±0.10 


25.36±0.09 


24.62±0.06 


23.52±0.27 


23.49±0.03 




22.50±0.03 


21 


31±0 


08* 


20.39±0.04 


20 


49±0 


18 


226 


25.40±0.11 


25.30±0.09 


25.48±0.11 


25.23±0.09 


25.17±0.08 


24.59±1.00 


24.75±0.08 


24.78±0.68 


24.37±0.11 


23 


41±0 


46* 


21.90±0.10 


21 


75±0 


36 


228 


<28.50* 


<27.70* 


<27.70* 


<26.80* 


26.17±0.16 




25.17±0.10 


25.17±0.67 


24.38±0.11 


23 


10±0 


40* 


21.14±0.06 


21 


23±0 


22 


234 


27.40±0.51* 


26.43±0.22* 


26.65±0.28* 


25.50±0.12* 


24.88±0.08* 


23.76±0.30 


23.52±0.08* 




22.76±0.03 


21 


38±0 


20* 


20.31±0.11* 


20 


57±0 


39 


236 


20.92±0.01 


20.73±0.01 


21.02±0.01 


20.59±0.01 


20.65±0.01 


20.32±0.02 




20.26±0.02 


19.94±0.01 


20 


05±0 


04* 


19.27±0.02 


19 


32±0 


04 


258 


24.01±0.04 


23.54±0.04 


23.70±0.03 


23.19±0.03 


22.68±0.02 


21.87±0.06 


21.85±0.01 


21.56±0.04 


21.11±0.01 


20 


61±0 


05* 


19.62±0.04 


19 


71±0 


04 


285 


24.60±0.18* 


24.65±0.15* 


24.63±0.24* 


23.88±0.14* 


23.70±0.12* 


23.32±0.28 


23.06±0.10* 




22.36±0.11* 


21 


80±0 


13* 


20.78±0.10* 


20 


84±0 


10 



NOTE. — Column description: (1) ID number of the object; (2) CFHT u* magnitude with its error; (3)- (4)- (5)- (6) Subaru Bj, g+ . Vj, r+ magnitudes with their errors; (7) CFHT i* magnitude with its error; (8) Subaru 
+ magnitude with its error; (9) HST/ACS F'814W magnitude with its error; (10) Subaru z+ magnitude with its error; (11) UKIRT J magnitude with its error; (12) CFHT K magnitude with its error; (13) NOAOKg with 
ts error. The values marked by * arc measured by our 3"-apcrturc photometry on the images. 
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TABLE 7 

COSMOS MULTI WAVELENGTH COUNTERPARTS OF THE SAMPLE 



ID 


FXJV N U V 


IRACl 


IRACl 


IRAC3 


IRACi 


MIPS 


1 




54.89±0.23* 


36.50±0.29* 


23.87±0.97* 


14.73±2.07* 


<0.15 


2 




41.44±0.19 


45.75±0.30 


28.07±0.99 


21.8±2.49 


<0.15 


3 




13.65±0.15 


23.62±0.25 


55.88±0.94 


166.89±1.98 


1.47±0.02 


4 




42.97±0.19 


48.78±0.29 


37.71±1.04 


29.34±2.19 


<0.08 


5 




66.44±0.19 


89.03±0.29 


91.57±0.99 


63.65±1.97 


0.87±0.02 


11 




25.34±0.17 


29.34±0.26 


25.07±1.00 


9.64±1.98 


<0.15 


13 




49.51±0.18 


51.39±0.26 


53.24±0.98 


75.84±2.04 


0.32±0.03 


16 




44.29±0.18 


32.54±0.28 


26.19±0.94 


18.39±2.26 


0.22±0.02 


18 




31.03±0.17 


24.94±0.28 


18.13±0.99 


14.03±2.15 


<0.61* 


20 




63.11±0.18 


44.73±0.27 


25.12±0.90 


17.46±2.06 


<0.08 


22 




57.27± 0.17 


64.90±0.29 


41.55±0.89 


32.84±2.16 


0.25±0.03* 


25 




49.17±0.19 


55.52±0.28 


38.82±0.99 


25.08±2.19 


0.18±0.05* 


26 




126.74±0.25 


108.17±0.29 


59.54±0.98 


44.75±1.92 


<0.15 


28 




32.46±0.16 


39.98±0.26 


46.73±0.91 


40.34±1.98 


0.15±0.02 


29 




9.81±0.16 


10.17±0.25 


6.50±1.00* 


<7.38* 


<0.15 


30 




74.26±0.19 


58.67±0.28 


33.32±0.91 


22.76±2.05 


<0.15 


31 


25.15±0.24 


47.86±0.16 


34.60±0.25 


24.98±0.81 


19.36±1.78 


<0.15 


32 




16.97±0.15 


32.41±0.24 


58.17±0.95 


98.81±1.81 


0.26±0.02 


34 




22.80±0.44* 


26.73±0.48* 


25.44±1.10* 


12.07±2.40* 


<0.26* 


36 




49.50±0.16 


43.80±0.27 


30.24±0.87 


21.89±1.97 


0.15±0.04* 


37 


23.55±0.15 


88.32±0.22 


123.39±0.32 


158.41±1.10 


249.52±2.18 


1.58±0.39 


38 




30.47±0.17 


28.00±0.28 


29.99±0.95 


24.94±2.35 


0.13±0.06* 


39 




74.19±0.21 


66.35±0.32 


34.74±0.96 


25.29±2.40 


<0.08 


52 


24.97±0.21 


59.84±0.21 


41.39±0.28 


36.37±1.09 


28.67±2.04 


0.16±0.02 


70 




17.58±0.16 


20.57±0.26 


29.83±1.04 


21.32±2.17 


0.13±0.03* 


202 




24.81±0.14 


22.02±0.24 


15.90±0.84 


6.46±1.87 


<0.15 


219 




98.71±0.23 


86.12±0.32 


48.50±1.02 


38.60±2.26 


0.12±0.01 


224 




48.09±0.17 


41.65±0.26 


25.56±0.88 


18.53±2.07 


<0.15 


226 




16.94±0.15 


30.93±0.28 


53.89±0.89 


104.92±2.16 


0.54±0.02 


228 




39.72±0.18 


47.82±0.27 


36.06±0.99 


23.46±2.03 


<0.15 


234 




50.19±0.19 


45.78±0.26 


33.95±1.03 


14.69±2.00 


<0.15 


236 




99.23±0.20 


146.61±0.36 


224.02±1.04 


324.18±2.50 


0.98±0.04 


258 


24.97±0.15 24.32±0.09 


84.11±0.23 


59.75±0.30 


37.90±1.06 


31.78±2.12 


<0.15 


285 


24.85±0.18 


34.06±0.17 


33.78±0.26 


33.68±0.95 


32.70±2.03 


0.13±0.03* 



NOTE. — Column description : (1) ID number of the object; (2)-(3) GALBX FUV and NUV magnitudes with their errors; (4)- (5)- (6)- (7) Spitzer/IRAC 4-channel (3.6, 
4.5, 5.8, and 8.0 /^m) fluxes with their errors; (8) Spitzcr/MIPS flux at 24/^m with its error. The values marked by * arc measured by our 3"-apcrture photometry on the 
images. 



